Hidden symmetry and knot solitons in a charged two-condensate Bose system 
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We show that a charged two-condensate Ginzburg-Landau model or equivalently a Gross- 
Pitaevskii functional for two charged Bose condensates, can be mapped onto a version of the nonlin- 
ear 0(3) (j-model. This implies in particular that such a system possesses a hidden 0(3) symmetry 
and allows for the formation of stable knotted solitons. 



For over forty years there has been a wide interest in 
condensed matter systems with several coexisting Bose 
condensates [0. Here we shall investigate the physi- 
cally important example of two charged condensates to- 
gether with their electromagnetic interaction. This sys- 
tem is described by a Ginzburg-Landau model with two 
flavors of Cooper pairs. Alternatively, it relates to a 
Gross-Pitaevskii functional with two charged condensates 
of tightly bound fcrmion pairs, or some other charged 
bosonic fields. Such theoretical models have a wide range 
of applications and have been previously considered in 
connection of two-band superconductivity. Indeed, these 
models describe superconductivity in transition metals 
The presence of two condensates has been observed 
in experiments on Nb, Ta and V as well as in iV&-doped 
SrTiOz More recently the renewed interest to two- 
gap superconductivity was sparked by discovery of the 
the two-band superconductor with surprisingly high crit- 
ical temperature MgB2 Q|. It has also been argued in 
lH that under certain conditions liquid metallic hydro- 
gen might allow for the coexistence of superconductivity 
with both electronic and protonic Cooper pairs. In a liq- 
uid metallic deuterium a deuteron superfluidity may sim- 
ilarly coexist with superconductivity of electronic Cooper 
pairs ^ (see also 0). 

Here we shall be particularly interested in an exact 
equivalence between the two-flavor Ginzburg-Landau- 
Gross-Pitaevskii (GLGP) model and a version of the 
nonlinear 0(3) cr-model introduced in We expose 
this equivalence by presenting an exact, explicit change 
of variables between the two models. The model in ||] 
is particularly interesting since it describes topological 
excitations in the form of stable, finite length knotted 
closed vortices The equivalence then implies that a 
system with two charged condensates similarly supports 
topologically nontrivial, knotted solitons. Previously it 
has been argued that these topological defects could play 
an important role in high energy physics ||l^. The 
purpose of the present Letter is to discuss the condensed 
matter counterparts. 



A system with two electromagnetically coupled, oppo- 
sitely charged Bose condensates can be described by a 
two-flavour (denoted by a = 1,2) Ginzburg-Landau or 
Gross-Pitaevskii (GLGP) functional. 
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where we take 
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Here we shall consider the general case where the two con- 
densates are characterized by different effective masses 
rria, different coherence lengths = S/ \/2maba and 
different concentrations Na —< >= ba/ca Q. 

The properties of the corresponding model with a sin- 
gle charged Bose field are well known. In that case the 
field degrees of freedom are the massive modulus of the 
single complex order parameter and a vector field which 
gains a mass due to the Meissner-Higgs effect. The im- 
portant property of the present GLGP model is that the 
two charged fields are not independent but there is a 
nontrival coupling which is mediated by the electromag- 
netic field. This implies that we have a nontrivial, hidden 
topology which becomes obscured when we represent the 
model in the variables ([|). In order to expose the topo- 
logical structure we introduce a new set variables, in- 
volving a massive field p which is related to the densities 
of the Cooper pairs and a three-component unit vector 
field n. The important feature of these new variables is 
their gauge invariance: Neither the relative phase of the 
condensates ^'i and nor the gauge field A enters in 
the free energy functional when represented in these new 
variables. 

We start by introducing variables p and xi,2 by 
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where the complex Xa = |Xa|e*'''° are chosen so that 
IXiP + 1X2 P = 1- The modulus p then has the following 
expression: 
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In terms of the variables p and Xa the free energy density 
in (111) reads as follows: 
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where we denote — |f ^fe, and d (without an index) is 
the ordinary V operator. The standard gauge invariant 
expression for the supercurrent density 

3 = — wid^i - 'j'la^n - -^28%} 
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becomes in these new variables 
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where 



j = i[xidx*i - Xidxi - X2dx*2 + X*2dx2]- 



We introduce a gauge-invariant vector field C, directly 
related to the supercurrent density by 



(5) 



We then rearrange the terms in as follows: We add 
and subtract from (||) a 
the following expression 



and subtract from (H) a term p^'i^ and observe that 
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is also gauge invariant. Indeed if we introduce the gauge 
invariant field 



n = (x , S x), 



(7) 



where x = (xii X2) and a are the Pauli matrices, then n 
is a unit vector Inl = 1 and we can write (0) as follows, 



fc2 2 

h p 



;2i 



\dxi? + \dx2?-' 



= \h'p\dnf (8) 



Consider now the remaining terms in (|^). For the mag- 
netic field we get 



B = rotA 
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where rot j can be written in terms of the unit vector n 
as follows: 



1 



rotj ~ -n • diXY x djii 



Combining these we then arrive at our main result: The 
GLGP free energy density becomes 



F^^{dnf + h^{dpY + 
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where we identify a version of the nonlinear 0(3) tj-model 
introduced in M 
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in interaction with a vector field C. With this we 
have completed the mapping between the two-condensate 
GLGP model (|l|) and (|l^, which is an extension of the 
model (|Tl|) introduced in We emphasize that this in- 
volves an exact change of variables between the GLGP 
model and (^0|). This change of variables in particular 
eliminates the gauge field and as a consequence the final 
result involves only the physically relevant field degrees of 
freedom that are present in the two-condensate system. 
Note in particular the appearance of the mass for 
which is a manifestation of the Meissner effect; 
the London magnetic field penetration length is = 
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We also emphasize 
that the contribution n • diii x 9,11 to the magnetic field 



term in (10), ( |ll| ) is a fundamentally important property 
of the two-condensate system which has no counterpart 
in a single condensate system. Indeed, it is exactly due to 
the presence of this term that the two-condensate system 
acquires properties which are qualitatively very different 
from those of a single-condensate system: This term de- 
scribes the magnetic field that becomes induced in the 
system due to a nontrivial electromagnetic interaction 
between the two condensates. 

The potential term V depends only on Cooper pair 
concentrations and masses, and it is a function of the 
modulus p and the 713 component of the vector n only. In 
particular, we can write the mass term in ( pO| ) explicitely 
as follows: 



V = A + Bns + Cnl 
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Where 
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A = p^[Acim\ + Ac2m\ — biuii — 627712] 
B = p^[%C2m\ — %cim\ ~ 52TO2 + biuii] 
C — 4p^[cim^ + 021712] 



(13) 



This mass term determines the energeticahy preferred 
ground state value for 773, which we denote by 77,3. Ex- 
phcitely, 
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Thus the ground state value of n is a circle specified by 
the condition 773 = 773 on the unit two-sphere. This yields 
a uniform, unperturbed ground state for the condensates. 
Note that the ground state value only depends on the 
concentrations Na and the masses 777q, of the Cooper 
pairs. 

However, here we are mostly interested in topological 
defects where the unit vector n locally deviates from this 
ground state value. This then corresponds to the for- 
mation of a local inhomogeneity in the densities of the 
Cooper pairs. But in order to ensure a finite energy, at 
a certain distance from the inhomogeneity the vector n 
should return to the circle defined by (14). 

Indeed, the equivalence between the models (0) and 
( p^ ) reveals a hidden topological structure in (|l|) which 
has a number of important physical consequences. In 
particular, it leads to the understanding of vortex types 
which are allowed by a two-condensate system: The 
model (^l|) is known to admit topological defects that 
have the shape of stable knotted solitons and are charac- 
terized by a nontrivial Hopf invariant (see also 
comment |l5|). The simplest nontrivial solution to the 
equations of motion that follow from ( pT| ) was proposed 
in . In terms of vector n it forms a toroidal vortex ring, 
twisted once around its core before joining the ends. The 
equivalence between the two-condensate model (^ and 
( [lO| ) then implies that the GLGP model with two charged 
condensates not only possesses a hidden 0(3) symmetry, 
but should also allow for the formation of such knotted 
vortex solitons. The stability of these topological defects 
has been confirmed in extensive numerical simulations 
1^- |l^ ] (for video animation of the numerical simula- 
tions of the knotted solitons see the www-address indi- 
cated in Ref. [Q). Indeed, the knotted vortex solitons 
6.17. in a two-condensate superconductor should consist of 
finite-length stable closed vortices, which carry a nontriv- 
ial helical configuration of the magnetic field. These finite 
length closed vortices have properties which are substan- 
tially different from those of a loop which is formed by 
an ordinary Abrikosov vortex. In particular: the present 
vortices are protected against shrinkage by the third term 
in (10) [second term in (pT|)]. 

Since the variable 773 acquires the preferable value ( p^ ) 
in the ground state, the knotted solitons have the follow- 
ing form: At spatial infinity the unit vector n assumes 



a value on the circle which is defined by the condition 
77-3 = 773. We denote this value of n at spatial infinity by 
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= n(x^oo)). (15) 



At the center of a knotted soliton the unit vector then 
reaches a value which corresponds to the opposite point 
on the unit sphere. In the case of a nontrivial knotted 
soliton the curve where the vector n coincides with this 
opposite point value in general forms a closed loop or 
a knot. This is the core of the knotted soliton and we 
denote it by [pf 
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At the core the densities of the condensates are charac- 
terized by the following non- vanishing values |\l/i(no)P = 
2ii7V2; |«'2(no)P = ^Ni. In between the core where 
n = no and the boundary where n n°° the unit vector 
n in general rotates in a manner which is determined by 
the ensuing Hopf invariant ||l9| . Consequently a knotted 
solitons can be characterized by the number of rotations 
which are performed by the component of n, which is 
perpendicular to the axis defined by the boundary condi- 
tions n°° and no when we move around the knotted soli- 
ton by covering it once in toroidal and once in poloidal 
directions over any surface which is located between the 
core n = no and the boundary n = n°°. In particular, in 
the region between the core and the boundary the mag- 
netic field acquires a contribution from n with a helical 
geometry, B'* cx n • 9^11 x djU. 

A schematic plot of a toroidal-shaped knot soliton is 
given on Fig. 1. 
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FIG. 1. A schematic plot illustrating helical geometry of 
the vector n in a simplest toroidal knot soliton (for simplicity 
we neglected anisotropy effects connected with nonzero mass 
of the na component of the vector n). The toroidal hollow 
surface is a surface situated between the knot core no and 
the knot boundary n°°. On this surface the vector n has a 
constant projection to the axis defined by the boundary con- 
ditions n°° and n. The spiral lines on this surface indicate the 
lines where the vector n has some particular constant position 
on the unit sphere (the picture schematically shows a toroidal 
knot soliton which, in terms of the vector n, was twisted twice 
around its core before joining its ends). If we proceed along 
an arbitrary path over this surface covering the knot soliton 
once in toroidal direction and then once in poloidal direction, 
the vector n will perform n and m rotations correspondingly, 
according to the Hopf invariant which characterizes the given 
soliton. In a general case the fluxtube can also be knotted 
(see also animations available at http address [9] ). 

Finally we comment on the length scales that are 
present in our system. For simplicity we assume that the 
two condensates have equal coherence lengths = ^2 = 
^. In this case one concludes immediately that (|l^) has 
two different length scales. These are A which is the mass 
of the field C, and ^ which is the mass of the component 
n3 and relates to the coherence lengths of condensates. 
In the typc-I limit which corresponds to A << ^, and in 
the type-II limit where A >> ^ the characteristic sizes of 
knotted solitons are then determined by different factors: 
In the type-I limit the size is determined by an interplay 
of the terms {dn)^ and {n-diiix djii)'^ which implies that 
the size is of order A In the type II limit the size 

of the soliton is determined by the large mass of and 
thus it is of order of \/Aj. 

In conclusion, we have argued that charged two- 
condensate Bose systems possess a number of interesting 
properties which are qualitatively different from those of 
a single condensate system. In particular they can sup- 
port stable knotted solitons as topological defects. In- 
deed, we find it remarkable that a GLGP functional for 
a two-component charged superfiuid becomes intimately 
related to the model introduced by one of the authors 
in which has been previously found to be relevant for 
strong interaction physics [p|jlC[| . This is a manifestation 
of the universal character of the model, which appears to 
describe a wide range of systems despite the differences in 
their physical origins. As a consequence the possibility 
of experimental investigation of the formation, proper- 
ties and interactions of knotted solitons in approproate 
superfluids ( e.g. Nb, Ta, V , iV&-doped SrTiOs) 

along with numerical simulations, could then shed light 
even on the properties of similar objects in nonabelian 
gauge theories of fundamental interactions. 
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